Dynamics of a Tonks-Girardeau gas released from a hard-wall trap 
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We study the expansion dynamics of a Tonks-Girardeau gas released from a hard wall trap. Using 
the Fermi-Bose map, the density profile is found analytically and shown to differ from that one of 
a classical gas in the microcanonical ensemble even at macroscopic level, for any observation time 
larger than a critical time. The relevant time scale arises as a consequence of fermionization. 
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The Tonks-Girardeau (TG) regime Q is that of ID 
impenetrable bosons, which is relevant to atom waveg- 
uide experiments with low densites and temperatures 
and large scattering lengths @. Under such conditions 
the radial degrees of freedom are reduced to zero-point 
oscillations, resulting a ID effective system, as it has 
been demonstrated in several experiments [j| ■ Some re- 
markable studies dealing with dynamics in this regime 
have shown the limits of mean-field theory in splitting- 
recombination processes 0, the spatial focusingof the 
probability density through Talbot oscillations [a], and 
existence of dark and grey solitons in a toroidal trap 
0. During the ID expansion of a harmonically confined 
Tonks gas, fermionization of the system was observed in 
the momentum distribution A parallel study point- 

ing out the "reciprocal" character of the Fermi-Bose du- 
ality, was recently carried out in for a fermionic TG 
gas which undergoes a dynamical bosonization. In the 
mean time, and motivated by the experimental build-up 
of square well 0] and hard wall optical box traps [ijj a 
growing interest has been developed concerning low di- 
mensional Bose gases trapped in such geometries ^3, El 
after the seminal paper by Gaudin Q . However, most of 
the studies have dealt with the gas within the trap and, 
at variance with the case of a harmonic confinement, only 
the single-particle evolution has been considered, in the 
field of ultracold neutron interferometry and diffrac- 
tion in time fl6| . 

In this Letter we account for a detailed study of the 
expansion dynamics of a TG gas, after switching off 
the confining hard wall potential and show the devia- 
tion from the associated classical gas in the microcanon- 
ical ensemble. In such a regime the Fermi-Bose (FB) 
map P, 0, 0, |2(J gives the many body wavefunc- 
tion of N strongly interacting bosons from the one of 
a free Fermi gas with all spins frozen in the same direc- 
tion. In order to do so, it suffices to apply the "anti- 
symmetric unit function" A = Tli<j<k<N s 9 n ( x k — Xj), 
as iPb(xi, ■ ■ ■ ,x N ) = A(x 1 , . . . ,x n ~)iIj f (xi, . . . ,x N ). The 
Fermi wavefunction, antisymmetric under permutation 
of particles, is built as a Slater determinant, with one 
particle in each eigenstate of the trap 



One more advantage of the general FB mapping is that 
it holds for time dependent processes (governed by one- 
body external potentials), since the A operator does not 
include time explicitly O, H3| , this is, Fermi statistics 
holds under time evolution. The glaring upshot is that 
as far as local coordinate distributions are concerned, to 
deal with the manybody Tonks gas it suffices to work out 
the single particle problem, since \iPb(xi, . . . , XN',t)\ = 
|^f(xi, . . . , xjv; t)\ 2 . In particular, from the involutiv- 
ity of the A operator (A 2 — 1) and the fact that 
((f>n\U U\(f> m ) — S nm , where U is the time-evolution op- 
erator, it follows that the time-dependent density profile 
can be calculated as 201 



p(x,t) = N \i/jb(x,X2, ■ • • 7 XN;t)\ 2 dx 2 ■ ■ -dx 
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Single eigenmode dynamics. Motivated by the Eq. , 
in this section we tackle the problem of studying the time 
evolution of the n-th eigenstate of a hard wall trap. As 
it is well-known they have the general form (f) n (x,t = 

0) = y^sin(n7rx/L)x[o,L] (x), with n E N where the 
characteristic function can be conveniently written as a 
difference of Heaviside functions, X[o.l]{x) = 0(L — x) — 
Q(—x). Concerning the dynamics we consider the free 
evolution under the propagator 
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after suddenly switching off the trap at zero time. Using 
the superposition principle and introducing p n — hnn/L, 



4> n (x,t) 
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FIG. 1: Density plot of the probability density for the dynam- 
ical evolution for the fifth eigenstate ( 8r Rb atom, L — 80/^m). 




100 200 300 400 

t(rns) 



500 



FIG. 2: Density plot of the probability density p for a Tonks 
gas composed of N = 5 atoms of 87 Rb (L — 80fim). 
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and w(z ) = e z er fc(- iz) the so called Faddeyeva func- 
tion |2lJ. After |22(, M(x, k, r) has been named the 
Moshinsky function. Similar expressions have been de- 
rived in the field of ultracold neutron interferometry , 
and discussed recently in the context of diffraction in time 
[lif . When the particle is trapped in the n-th mode, and 
at short times after being released, the probability den- 
sity presents n maxima, but the central ones tend to fade 
away with time as shown in Fig. ^ Indeed, the general 
structure of the n-th eigenstate (n > 1) under evolu- 
tion presents a bifurcation in two main branches after 
the semiclassical time t n — mL/(2p n ) = mL 2 / (2nirh) . 

Tonks gas dynamics. Remarkably, as stated above, the 
calculation of the density profile is possible from Eqs. 

and ©. As illustrated in Fig. [3 the density 
profile exhibits at short times a peaked structure with 
the number of maxima equal to that of particles. 

This spatial antibunching results from the underly- 
ing fermionization characteristic of the Tonks regime. 
Indeed, the two-particle local correlation, gi = 
(* t (a;) 2 *(a;) 2 )/(*I' t (a;)*(2;)) 2 , was shown to tend to zero 
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FIG. 3: a) Dependence of the variance on the number of 
87 Rb atoms for the density profile of TG gas confined in a 
trap, b) Time dependence of the variance of the density profile 
for different number of 87 Rb atoms. Filled circles mark the 
corresponding tjy and L — 80/xm. 



even for non-zero temperatures |25| . However, the visibil- 
ity of the interference pattern is lost both with increasing 
time and the number of bosons under consideration. 

A measure of the roughness of the density profile is the 
root mean square of the density profile weighted with 

itself, a 0r , 
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N/PN 



, with pn the density 
profile normalised to one particle. In fact, at time equals 
zero and as a function of the particle number this mea- 
sure decreases monotonically as shown in Fig. [3^, in 
agreement with the resolution of the identity within the 



box, Zirnjv- 



dP — 1[0,L] 



see Eq. ©. The time evolution 



is richer in structure, presenting an initial positive slope, 
see Fig. |3Jd, to finally start to decrease after t ~ t^, 
reaching zero asymptotically. 

After a transient regime, where the structure of the 
peaks varies in time, the roughness of the density pro- 
file becomes monotonically decreasing with both time 
and number of particles. In particular, the fastest 
components will be associated with the highest ex- 
cited state (the Fermi level in the dual system of non- 
interacting spin- frozen fermions), and therefore with 
quasi-momentum given by ±pjy = ±hNir/L. Actually, 
these components govern the width of the expanding 
cloud for t ^ tjsf. Figure 0] shows the variation in time 
of the full width at half maximum (FWHM) for differ- 
ent number of particles. The upshot is that deviations 
from a linear dependence on time are observed right after 
switching off the trap, similarly to the results reported in 
8] for the harmonically confined Tonks gas. Such devia- 
tions disappear with increasing number of particles and 
already for TV ~ 100 are bellow the millisecond time scale. 
From its definition, one can conclude that tjy, the time 
necessary for a classical particle to leave the trap when it 
is initially located in the center of the box and moves with 
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FIG. 4: Time dependence of the FWHM of a cloud of 1 
atoms released from a trap 80/im in width. 
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FIG. 5: Measure Dm for a fixed density, do = 1 atom/ fim 
and increasing number of particles ( 87 Rb). Filled circles mark 
tN and verify Djv(tiv) — 0.08. 



the momentum ppj , arises as a consequence of fermioniza- 
tion. Nevertheless, it is clear from the transient features 
of the single-particle solution that a self-similar expan- 
sion does not occur for a TG gas released from a hard 
wall trap. This fact contrasts with the harmonic case, 
pointing out the relevance of the confining geometry. 

Comparison with a classical gas. Next we consider 
a classical gas of distinguishable particles interacting 
through a contact potential. Since for the TG gas the 
total energy (E = J2^=i E „,E n = h 2 n 2 it 2 / {2mL 2 )) is 
fixed once the number of particles N in the initial trap is 
specified ^jj , it is natural to deal with a microcanonical 
ensemble in the classical limit. This leads to a distribu- 
tion in phase space of the form: 
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Y[Q{ Xi )Q{L- Xi ), 
(6) 



reads 



where the normalization 
NT(N/2)/(2rmr N / 2 L N g N - 1 ), 
function and g = \/2mE. 
orem, momentum 
change), and the symmetry 
of particles in the microcanonical 



A/o 

r(z) is the Gamma 
Using Liouville's the- 
conservation (but for inter- 
under permutation 
ensemble, one 



finds p cl (x,p,t;E) = p c /(x - pt/m,p,0;E), where 
x = {xi, - ■ ■ ,x N }, and p = {p%, ■ ■ ■ ,Pn}- To 
obtain the density profile, an integration over 
{x2, ■ ■ ■ , Xpf,px, . . . ,pjv} variables has to be performed. 
Note that each of the integrals over Xi contributes 
exactly L. Moreover, all the integrals over momenta can 
be carried out by using the generalisation of spherical 
polar coordinates to a hypersphere in a ./V-dimensional 
space subjected to the constraint YliPi = Q 2 = ^ m E 
[25.l2^| . which fixes the radius g of the hypershell. Since 
it is always possible to choose one coordinate of the form 
p = Pi = Q cos £, it follows that 

p c l(x,t;E) = M ( d£sm N ~ 2 £Q(x-tgcos£/m) 
Jo 

x6(L + tgcos^/m — x), (7) 
where TV = NT(N/2)/{^LT[(N - l)/2]) is the normal- 



isation constant. 

For a systematic comparison between the Tonks and 
microcanonical gas, we introduce the time-dependent 
measure Djy(t) = J dx\p{x,t) — p c i(x,t)\/N , which is 
plotted in Fig. [5] for different number of particles. As a 
result, a clear time scale given by ijv can be established, 
such that for small times (t/tfq < 1) both quantum and 
classical profiles essentially coincide (D^ < 0.1). 

In particular, the microcanonical model reproduces for 
short times and low number of particles the quantum 
profile in a coarse-grained fashion, neglecting inteference 
(see FigEJi). It is precisely for low TV when the differences 
between the initial quantum spatial distribution and the 
uniform classical one, are greater. For larger number of 
particles the quantum profile tends to "resolve the iden- 
tity" within the trap, reaching a uniform constant distri- 
bution when N — > oo. This is, the quantum result is a 
staggeringly flat distribution but for the smoothed edges, 
and can be understood as a consequence of the discrete 
spectrum of the free Hamiltonian confined in the sub- 
space [0,L], with equally spaced quasimomenta p n (See 
FigHD- Indeed, in the limit t/t N < 1 and N > 1, the 
profile is well described by the characteristic function of 
the form Q(x + pFt/m)<d(L + ppt/m — x)/(L + 2ppt/m) 
(Fermi- hat), which expands with the momentum of the 
Fermi level in the dual system, p^. Similar profiles have 
already been reported for an effective ID weakly inter- 
acting fermionic gas confined in a harmonic trap [2r| . 

By contrast, it can be proved that the momentum dis- 
tribution of the microcanonical gas is of the form [24| 
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which in the reduced variable p/g, for large AT, tends 
to be a Gaussian with zero mean and 1/N variance. It 
follows that the density profile becomes also Gaussian in 
this limit, in clear disagreement with the TG profile, as 
shown in Fig. [S]for t > t^. 

In the thermodynamic limit, achieved at a given ob- 
servation time and fixed density do = N/L by taking 
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FIG. 6: Microcanonical (dashed lines) and quantum (contin- 
uous line) density profiles 100ms after expansion for: a) fixed 
length box L = 80/xm, and N(t N )=W (137ms), 25 (55ms), 50 
(27ms), and 100 (1ms), from top to bottom, b) fixed density 
d = latom//im, and N(t N ) =100 (21ms), 250 (54ms), 500 
(107ms), and 2000 (428ms), from bottom to top. 

N — > do, both the quantum and microcanonical density 
profiles become indistinguishable of one another. Since 
for fixed density tjy = mN / (2Trhd^) , the agreement is 
reached more easily at short times after switching off the 
trap. Figure shows that at an observation time of 
100/is, 2000 particles suffice. Deviation from the classi- 
cal gas increases with time which acts as a lense pointing 
out the difference in the underlying momentum distribu- 
tions. Indeed, taking note of the experimental number of 
atoms in generally between 500-3500, a rectangular 
profile is to be expected in the TG regime both in the 
trap and after switching it off, being overall described by 
the Fermi-hat model. 

Discussion. Taking advantage of the Fermi-Bose map, 
the dynamics of a strongly interacting many-body system 
can be worked out once the time evolution of the single- 
particle eigenstates of the confining hard-wall trap is 
known. In this way, we have accounted for the first st udy 
of a many-body Moshinsky shutter problem [22], |2J], |2J| 
including interactions. Right after switching off the trap, 
the profile evolves through a transient regime where the 
overall behaviour is analogous to the one of a classical 
gas in the microcanonical ensemble. 

However, at times t > t/v — ml? / (2Nirfr) and in the 
same conditions (N,L), a regime with uniform density 
is observed in the Tonks gas, clearly deviating from the 
classical bell-shaped profile. The origin of this main fea- 
ture can be traced back to the discrete spectrum of the 
original trap potential. Finally, it is noteworthy that the 
results obtained in this Letter equally hold for a nonin- 
teracting Fermi gas due to the involutivity of the anti- 
symmetric unit function which entails the ID Fermi-Bose 
duality. 
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